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Abstract 

We introduce an extension of cone splines and box splines to fractional and complex orders. These 
new families of multivariate splines are defined in the Fourier domain along certain s-dimensional meshes 
and include as special cases the three-directional box splines 0 and hex splines fTH] previously considered 
by Condat, Van De Ville et al. These cone and hex splines of fractional and complex order generalize the 
univariate fractional and complex B-splines defined in [n m and investigated in, e.g., mm ■ Explicit 
time domain representations are derived for these splines on 3-directional meshes. We present some 
properties of these two multivariate spline families such as recurrence, decay and refinement. Finally it 
is shown that a bivariate hex spline and its integer lattice translates form a Riesz basis of its linear span. 

Keywords and Phrases: Cone splines, box splines, s-dimensional mesh, hex splines, (fractional) 
difference operator, fractional and complex B-splines 

AMS Subject Classification (2010): 65D07, 41A15, 42A38 


1 Introduction 


We introduce an extension of cone and box splines to fractional and complex orders. Given a direction matrix 
or knot set, both types of multivariate splines are defined in the Fourier domain as ridges along the directions. 
These new families include as special cases the three-directional box splines introduced by Condat and Van 
De Ville [3] as well as the hex splines defined by Van De Ville et al. m • Whereas the latter multivariate 
splines provide integer smoothness along a direction, the new families allow for fractional smoothness along 
different directions thus defining continuous families (with respect to smoothness) of functions. In the case 
of complex orders, cone and hex splines have an additional phase factor automatically built in that allows - 
just as in the one-dimensional setting - the enhancement of frequencies along the directions defined by the 
knot set. The importance of phase has been well known since the experiment by Oppenheimer and Lim m 
For this reason, hex splines of complex orders may be potentially employed for multivariate signal or image 
analysis. 

This paper, however, focuses on the mathematical foundations of cone and hex splines of fractional and 
complex orders. In the current setting, these two spline families generalize the univariate concepts of frac¬ 
tional and complex B-splines as defined by Blu and Unser respectively, Forster et al. JB| to higher 
dimensions. Previously, an extension of univariate complex B-splines to higher dimensions was proposed in 
m- This approach is based on a generalized version of the Hermite-Genocchi formula and produces the 
analog of multivariate simplex splines (of complex orders). Our approach differs also from that undertaken 
in [8j where the multivariate extension was done via ridge functions. The focus of our setting lies more on 
the geometry of the knot set. 

The structure of this paper is as follows. In Section [2j we briefly review the definitions of cone and box 
splines and present those properties that will be relevant in subsequent sections. The next section deals with 
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cone splines on 3-directional meshes. There, we extend the results obtained in [3] and obtain an explicit time 
domain representation of such cone splines on 3-directional meshes with different weights along the various 
directions. Hex splines with different weights along the given directions are defined in Section [3] and it is 
shown that the hex splines in mm are again special cases. Fractional cone and hex splines are introduced in 
Section [5] While they are defined in the frequency domain, explicit formulas for them in the time domain are 
presented. We show that fractional cone splines in R s constructed from s knots obey a recurrence formula 
that is similar to that derived by Micchelli m for polynomial cone splines. In the final section, we derive 
decay rates for hex splines and also show that hex splines are refinable functions. We also show that with 
modest assumptions placed on the knots, the set V consisting of a bivariate hex spline and its integer lattice 
translates form a Riesz basis for span(V). We conclude the paper by extending the definitions of fractional 
splines to introduce cone and hex splines of complex order. 


2 Review of Cone and Box Splines 


The R-spline of Curry and Schoenberg |4j is defined as follows: 

N 0 , n {x) := n [x 0 , ■■■, x n \ (• - x)+ _1 (1) 

where xo < • • • < x n with x$ < x n are the knots, x + = max(0, x) is the truncated power function, 
and [xo,...,x„] is the divided difference operator. Multiplication by n normalizes the spline so that 
S N 0 , n (x) dx = 1. 


2.1 Cone Splines 


It is natural to define multivariate splines by first generalizing the truncated power function that appears in 
0. We have the following definition: 

Definition 2.1. Let X = {x 1 ,.--) 3 ^} <z R s be a linearly independent set of vectors with cone(X) = 
{fix 1 + • • • + t s x s | fi,..., t s ^0} and M = [x 1 • • • x s ]. Then the cone spline C(x | X) is defined as 



|det(M)|-\ 

0 , 


x G cone(X); 
otherwise , 


( 2 ) 


and inductively for additional nonzero knots by 


C (x | X u 


rCO 

{x s+1 }) = C( 
Jo 


x — tx' 


.8 + 1 


X) d t. 


Definition 
span(X) = 


2.1 


leads to several properties (see [2]) of cone splines. For X = {x 1 ,... , x™} <z R s \{0}, 
. The first four properties appear in [2] and the last can be found in [5]: 


n ^ s. 


• supp(C(x | X)) = cone(X) := (fix 1 + • • • + t n x n | ti ,..., t n ^ 0}. 

• C(x | X) > 0 for x G cone(X). 

• C(x | X) is a piecewise polynomial function of total degree n — s. 

• Let p be the minimum number of knots that can be removed from X so that the resulting set does not 
span R s . Then 

C(x | X) G C p ~ 2 (IT). (3) 
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• For all continuous functions g : R s —> ffi. with compact support, 



C(x | X)g(x) dx = 



(4) 


As the set Co(R s ) of continuous functions with compact support contains the set C® of infinitely differentiable 
functions as a dense (in L p , p e [1, oo)) subset and the latter is dense in the Schwartz space <S>(M S ), we can 
interpret the cone spline C{» | X) as a tempered distribution by setting 


<C(.|X),«?>:= f 5 (i> 
J[0,co) Tl Vfc=l 


dt, \/geS(R s ). 


We can use (|2.l|) to see that the Fourier transformation of the cone spline, considered as an element of 

1 


<S'(M S ), is given by 


n~- x) | 


My + 71-5 (<W,X fc >) 


(5) 


where f([ represents the direct or tensor product of distributions. Indeed, for any g e <S(R S ), 

/c(u> | X), cA = <C(x | X), = f gl Vi fc x fc ) dt = <H, 5 (Xt)> = /H(X T W ), 5 

\ / J[0,oo)» Vfc=l J ^ 


=/n^«x fc ^», 5 )=( n 


\k= 1 


\k =1 


i<w,x fe > 


+ nS «w,x fe » ) ,g 


where we set Xt := 2 and where H and H denotes the one-, respectively, n-dimensional Heaviside 

k= 1 

function. 

Remark 2.1. In the sequel, we need to consider direct products of generalized functions S andT of the form 
S m and S m T n , m,n e N. These are defined as follows JWj . 

(S m , 0(x)> := <S Xl , <5 Xa , ■ ■ ■ , cfix ,,..., x m )) ■■■)), 

where the subscript on S indicates on which variable S acts. Similarly, one defines 


<■ S m T n , 0(x, y)> := (S m ,(T n , 0(x, y» , 

where S m , respectively, T n acts on the m variables x = (ay,..., x m ), respectively, the n variables y = 


The distributional relation § can also be used to show that if A is any nonsingular s x s matrix, then 

C (x | (Ax 1 ,..., Ax"}) = (A -1 x | X) . (6) 

Finally, Micchelli m proved the following recurrence formula obeyed by cone splines: 

1 " 

C’(x | X) = ——— ^ AfcC (x | X\ {x fc }) (7) 

n s k= i 

n 

where x = 2 AfcX fc . Two cone splines are plotted in Figure 1 

k= 1 
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Figure 1: Two cone splines and their support. Let x 1 = (1,2) T and x 2 = (1,— 1) T . Then the knot set is 
X = lx 1 ,x 1 , x 2 , x 2 ,x 3 } and the top spline is a piecewise cubic polynomial. The spline on the bottom is 

piecewise linear with knots x 1 = l) T , x 2 = l) T , and x 3 = (1,0) T . 


2.2 Box Splines 


Multivariate box splines are defined in a manner similar to that for cone splines: 


Definition 2.2. Let X = {x 1 ,..., x s } c K s be a linearly independent set of vectors. Let M = [x 1 • • • x s ] 
and denote by Z(X) the zonotope {fix 1 + • • • + f s x s 1 ti,..., t s e [0,1]}. Then the box spline .B(x|X) is 
defined as 


f |det(M)|-\ xeZ(X); 
[ 0, otherwise, 


and inductively for additional nonzero knots by 


B (x | X u 


,*+1 


})= fs(> 

Jo 


tx 


8+1 


X) dt. 


Definition (2.2) gives rise to several properties (see [2]) of box splines. For X = {x 1 
n ^ s, span(X) = R s (see [2j): 


s \{0}, 


• supp(B(x | X)) = Z(X) = {fix 1 H-h t n x n | ti,... ,t n e [0,1]}. 

• B(x|X) > 0 for xe Z(X). 
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• B(x | X) is a piecewise polynomial function of total degree n — s. 

• For all continuous functions g : R s —> R. with compact support, 



S(x | X)g(x) dx = 



(8) 


As B(» | X) e L 1 (R S ), we interpret B(» | X) as a regular tempered distribution acting on Schwartz functions 
via Q. Employing similar arguments to those given in the derivation of the Fourier transform of a cone 
spline, replacing the n-dimensional Heaviside function by the characteristic function of the unit n-cube, one 
shows that B{» | X) is given by 

(l - e - i <"’ xh >j 


b(«|x)= n 

Writing the above equation in the form 


k=l 


i<w,x fe > 


(9) 


B(u\X)= H (l 


k= 1 


3 -*<u,,x fc >\ . 


Hi 

k=l 


Jj «<w,x fe > 


-no 

k =1 


7b / i 

3 -i<u,,x fc >\ . FT [ _£_ 


+ ttS «w,x fe » 


= : VxC'(wlX), 


( 10 ) 


where we set 


n 

V x := (l - e _i <“’ xfc >) 


k =1 


and used the fact that ^1 — e 1< A' X >j . $ (<(o;,x fc )) = 0 (in 5'(K S )). Taking the inverse Fourier transform of 

(101 and noting that Vx is the Fourier transform of the multivariate backwards difference operator pAy Vx 
defined by 

V x := [I Vx ’ 

xeX 

where the backwards difference operator in the direction x for a function / : R s —> M is given as 


V x / :=/—'/(• — x), 


we obtain the well known identity 

H(x|X) = -i T V x C(x|X). (11) 

r(n) 

We end this section with examples of box splines plotted in Figure [2j 


3 Cone Splines on 3-Directional Meshes 

A knot set X = {x^x^x 3 } cz ]R 2 will be called a 3 -directional mesh if x 3 ,x 2 are linearly independent and 
x 3 = x 1 + x 2 . 

We are interested in developing explicit formulas for cone splines on 3-directional meshes. In particular, we 
seek formulas for cone splines of the form 


C n (x | X) := C (ni , n2ins) ( x | X) := C x 


x 1 ,... x 1 , x 2 ,..., x 2 , x 3 ,... x 3 


n i times times n 3 times 
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Figure 2: Two box splines and their support. The top spline is a piecewise linear polynomial. The knots 
= (3,1) T and x 2 = (—2, — 2) T and x 3 = (—1,1) T . The spline on the bottom is piecewise cubic. If 


are 


= (|, l) , x 2 = —l) and x 3 = (—1, 0) T , the knot set is X = {3 




where rik e Z + := {n e Z I n > 0}, k = 1, 2, 3. 


For the special case n = (n,n,n), n e Z + and knot set X = { 


— fxbx 2 , x 3 


! } where 


x 1 = 


1 , 

VI j > 

2 


x 2 = 


3 , 

VI ) > 

2 


x 3 = X 1 + x 2 = ' 1 


Condat and Van De Ville [3| were able to derive a closed formula for C n (x | X): 

n t I Y -> 2 yi/n— l + k\ 1 

c n (x | x) = —j= 2_j 


\/3 —j V fc / (n - 1 - ft)!(2n - 1 + fc)! 


2|cr 2 | 

V3 


T— 1 — /c 


Xl - 


F 2 I 

V3 


2 m-l+fc 


( 12 ) 


It is straightforward to generalize (12) for the case when the knots of X are not repeated the same amount 


of times. In order to derive such a formula for n = (ni, 712 ,^ 3 ), we first need the following result for 
s-dimensional cone splines built from s knots. 


Proposition 3.1. Let X = {x 1 ,... ,x®]- c R s be a linearly independent set of vectors and suppose n = 


6 








(ni,..., n s ) e Zy. Then 


Cn(x | X) = |det(M)|- 1 [ [ 


<x,x fe y 


k\n k -1 


Ji 


(13) 


k =1 


where M = [x 1 • • • x s ] and {x 1 ,..., x s { zs t/ie dual basis to X. 

Proof. Let E = {e 1 ,..., e s } be the canonical basis for M s . Using ([5]), we have 

1 


C n (x|E) = [j 


fc = i \Ku,e k ) 


) n k S / , \Uk 

= n + ns • 


Applying the inverse Fourier transform to the above equation, noting that each factor in the n/c-fold direct 
product produces the truncated power function (xk)° + , we obtain the n/c-fold convolution product ( Xk)° + * 
• • • * (x k ) + - It is easy to establish that this product is ( Xk )" fc_1 . Thus we have 


C'n(x|E)= f] 


(Xk) 


n k -l 


i 1 ! (»*-!)!' 

Using ([ 6 ]), with A = M, and the fact that C n (x | X) = C n (x | {Me 1 ,..., Me s }), we can write 

1 


& W x >-P5(mUUm-'*|e)_ 


i n (xy 

|det(M)| (n k - 1)! 


and the proof is complete. 


□ 


We can use Proposition |3.1| with s = 2 to find a closed formula for cone splines whose knot set is a 3-directional 
mesh. We need the following lemma: 


Lemma 3.1. For n,m,£ e Z + and a;, a, 6 e R. with b < a < x, let 

CL,b) 

Then 


1 


(n— l)!(m—1)!(^— 1)! 


dit. 


n— 1 / 

k =0 ' 


m+k—l 
k 


f(u- 6 ) ri - 1 (u-a) m - 1 (x-w)^ 

J a 

' (a-b^-^x-a)™^- 1 . 


(n—l—k)\(m+£+k—l)\ 


Proof. The proof is a straightforward application of the result 


1 


T(a) 


J ( u-a ) /J \u-by 1 (x-u) a 1 du = (a — b) 1 1 ^ ( x-a) a+f3 t 2 F 1 ^1 - 7 ,/?;a + /3;-^ 77 ^ 

1 U (U) 
that appears in m p. 41] where b < a < x, a > 0, Re(/3) > 0, 7 e C. In our case, we insert positive integers 
7 = n, f3 = to, and a = l to write 


F(n,m,t) (Xj U, b ) 


1 

(n — 1 )!(to + £—!)! 


(a — b) n 1 (x — a) m+l 1 2 Fi ( 1 — n, m; m + l. 


(1 - n)k{m) k (~ 1 ) A 
(n — l)!(m + i — 1 )! (m + l)k k\ 


1 00 

m. A- f — 1 V 


(x — a) 


m+£+k— 1 


(a-b) 


a — b 


n—l—k 


(15) 
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The rising Pochhammer symbols in the 2 -Pi hypergeometric function are given by 

sk (n-l)\ 


(1 - n)k = (1 - n)(2 - n) ■ ■■ (k - n) = (-1) 
(m)fc = m(m +!)••• (m + k — 1) = 


(n—l — k)\ 
(m + k — 1)! 


(m — 1)! 


(m + £)k = 


(to + £ + fc — 1)! 
to + t — 1 )! 


Inserting these identities into (15), simplifying and noting [x — a) = (x — a) + gives the desired result. □ 
We are now ready to state a closed formula for cone splines on 3-directional meshes: 

Proposition 3.2. Suppose X c I 2 is the 3-directional knot set X = {x 1 ,x 2 ,x 3 } with x^x 2 linearly 
independent and x 3 = x 1 + x 2 = (1,0) T and, without loss of generality, x\ > 0. Then the cone spline 
C n (x I X) is given by 


Cn(x | X) = — < 


ni —1 


Y (ri2 + k— 1 \ _ 

Zj \ k ) (m-l-k 




fc =0 
n 2 — l 


n2+n3 + fc—1 


' l 2 a , / \ n 2 —l—k / 1 \ 

v /m + fc—1\_1_ (_I2 ) _ X 1 -l \ 

^ v fc / (n 2 — l—k)\(ni+n 2 + k — 1)! y y y 1 x 2 ^) 

Proof. From ([ 5 ]), we know that 

| X) = d (nun2) (w | {x 1 , x 2 }) • (— + 7 tS (wi)^ 

\iuii ) 


n 1 + ri3+k—l 


X 2 > 0 


X 2 < 0 


so that 


(ah') 713-1 

C„(x|X) = C (nijTl2) (x | {x 1 , x 2 }) * ^ Y. ^02) 


(n 3 - 1)! 

= Udet(M)| Y\ (nk - 1)M 

rx 1 

x <(u, x 2 ), x 1 )" 1- ^^, x 2 ) , x 2 )" 2-1 Oi - it)” 3 ' 1 du 

J—GO 


(16) 


Since x 1 + x 2 = (1,0) T , we know x 2 = (l — x\, —x 2 ) and from here we can easily show that the dual basis 

{i ‘* 2H {(40’(-k)l _ (17) 

and |det(M)| = x\. Inserting these identities into (16) and setting I\ = Hfc=i ( n k ~ 1)!^ gives 


C, 


,(x|X) = if J 


X1 ' (x\-l 

u — 1 


%2 


n i — l 


1 \ "2 — 1 


U - rX 2 {x\ — U ) 


"3 — 1 


dit. 


(18) 


If X '2 > 0, we have - L - s —x 2 < -^x 2 < Xi and (18) becomes 

















and we can use Lemma |3.1| to complete the proof. The case where X 2 < 0 is similar and we can use the 
continuity of the cone spline ^ for the case X 2 = 0 . 

□ 


Remark 3.1. Using ([ 6 ]), we can construct a formula for the cone spline with the general 3- directional mesh 


x 1 , x 2 , x 1 + x 2 . 


Remark 3.2. Setting n = (n, n,n) in Proposition 3.2 gives the cone spline (12) of Condat and Van De 
Ville. 


We use Proposition |3.2| to plot the cone splines shown in Figure [3] 



Figure 3: For X = {x^x^x 3 } where x 1 = , x 2 = and x 3 = (1,0) T , the cone splines 

plotted clockwise from the top left are C( lili3 )(x | X), C( 2 ,i.i)( x | X), C( 2j3i i)(x | X), and C'( 2 , 2 , 2 )( x | X). 


4 Hex Splines 


In this section we give a general formula for the so-called hex splines introduced by Van De Ville et. al. m- 
These splines were developed for applications that involve hexagonally sampled data. The hex splines can 
be defined by using box splines with the knot set X - = {x^yx 2 ,— x 3 } <z R 2 where x 3 = x 1 + x 2 . For this 
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knot set and n = (ri \. ri 2 - n-j) e Z 3 , our hex spline takes the form 



Figure 4: Support for hex splines (left to right) -B(i,i,i)(x | X ), -B( 2 ,i,i)(x | X ), -B(i, 2 , 2 )( x | X ), and 
B( 2 , 2 ,- 2 )( x I X - ) using the knots x 1 = , x 2 = , and —x 3 = —x 1 — x 2 = (—1,0) T . 

of the box spline: 


B n ( w|X“) = 


(1 - e -i<«. xi >j ni ^1 _ e -K^,x 2 >j " 2 ^ _ e ««,* 3 >y 

(^,xi» ni (z(u>, x 2 ))" 2 «u,,-x 3 >p 


^1 _ e -i<w, xl >^ ni ^ 


0 — 1( UUX 


)” 2 (e^>f) 




«U,, X l>n K ^ X 2 ))" 2 (,< Wl x 3»" 3 


1 - e-*"’* 1 )) \ ni / (l - \" 2 ( (e l<u; ' x3 > - l) ' 


(19) 


( 20 ) 


i^o^x 1 ) ) \ i(u;,x 2 ) J l i(uj,x 3 ) 

Using arguments from [T01 . Section 3.11] and [TTJ Theorem 7.1.16], we can rewrite each of the factors above 


as 


(l - e -*<“> x *>)' 
*<w,x fc ) 




for k = 1, 2 (and a slight modification when A: = 3) so that (20) becomes 

B n (o) | X") = (l - e- i < u ’- xJ >y 1 ' (l - e-^’X 2 ))" 2 (e z <“' x3 > - l)”' Cn(u: | X). 

Taking inverse Fourier transforms gives 

B n (x|X") = /(x)*C7 n (x|X) 
where /(x) is the inverse Fourier transform of 

/(«) = (l - e -*<“. xl >) ni (l - e -*<“.x 2 )^" 2 ^ e i<a;.x 3 > _ ^ " 3 


( 21 ) 

( 22 ) 

(23) 
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As noted in [3], if we can write (231 as 


/M = X c ki,k 2 e 


k> 


keZ 2 


where M = [x 1 x 2 ], then 


/( x ) = XI c fcl ,fc 2 5(x-Afk) 


keZ 2 


so that (221 becomes 


fl n (x|X-)= X c fcllfe C n (x-Mk). 


(24) 


keZ 2 


In the case where n = (ri, n, n), n e Z + , and x 3 = x 1 + x 2 = (1, 0) T , Condat and Van De Ville [3] found the 
following formula for the coefficients in (24): 


min(n+/ci ,n+&2 ,n) 

c ki,k 2 = X ( _1 ) 

k=max(ki ,k 2,0) 


k 2 +k 2 + k ( n 


n \ 1 n 


k — k\J V k — k 2 J \kj' 


(25) 


and combining (25) and (24) with the cone spline (121, they found the following formula for the hex spline: 

n n —1 


B„(x|X)= X C ^ 2 X 

ki,k 2 =—n k =0 


n — 1 + k 
k 


2x 2 

73 


+ ki - k 2 


I n—l—k 


Xl ~ 


(2 n — 1 + k)\(n — 1 — &;)! 
ki + k 2 | x 2 ki — k 2 


V3 + 2 


2 n—l+k 


(26) 


Figure [5] shows the hex splines of Condat and Van De Ville for n = 1,2. 
It is straightforward to generalize ([25]) to general n = ( 711 , 712 , 713 ): 


Proposition 4.1. For X = ■[x 1 ,x 2 ,x 3 } where x 3 = x 1 + x 2 = (1,0) T , M = [x 3 x 2 ] and n = ( 711 , 712 , 713 ) e 
Z 3 [231 can be expressed as 


f M = X c fcl , fc 2 e-*<“’ Mk > 


where 


n3+min(fei ,/c2,0) 

C ki,k 2 = X 

m=max(n3— n\ + fci ,713— n 2 + k 2 ,0) 


keZ 2 


Til 


712 


7ii — 713 — k\ + mj \n 2 — 713 — k 2 + mj \m 


n 3 


(- 1 ) 


fcl+^ 2+^1 


Proof. We follow the argument used by Condat and Van De Ville [3j. For the sake of notation, set z\ = 
e—*<“-x ) and 22 = x C Then [23| becomes 

/(«) = (1 - ^r 1 )" 1 (1 - z,:T 2 (*1*2 - ir 


Expanding each factor gives 


n 1 n 2 ri 3 

/m - 2 S S 

j=0 k =0 m=0 




n 3 -m-j 

1 


n 3 —m—k 

z 2 


Making the substitution k\ = j + m — n 3 and k 2 = k + m — n 3 and simplifying yields 


/M 


7T-3 n\+m—n 3 n2+'m—n 3 

XXX 

m —0 ki=m—n 3 k2=m—n 3 


n 1 


ni — n 3 + m — k± J \n 2 — n 3 + m — k 2 


712 



(-l) fcl 


+ fc 2 +m —k 1 




— &2 
2 
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Figure 5: (7 n (x | X) generator for B n (x | X) for n = 1 (top row) and n = 2 (bottom row). The functions in 
the top row are piecewise linear and the functions on the bottom row are piecewise quartic. The support 
sets for these splines are plotted (first and last) in Figure [4j 


We can interchange the order of summations by changing the limits on the inner two sums to k\ ,k 2 B Z 
and restricting the limits on m so that the binomial coefficient factors are not zero. These changes give the 
desired result for the Cfc 1; fc 2 . Finally, note that 


^ —fci ^-fc 2 _ -i(<dj,x 1 >fci+<i*>,x 2 >fc 2 ) _ -i(M T u)) T -k _ — i, 


■M k _ g—i<o 


i,Mk> 


(27) 

□ 


We can use Propositions 3.2 and 4.1 in conjunction with (24) to find explicit representations of hex splines. 
The hex splines corresponding to the second and third knot sets from Figure [4] are plotted in Figure [6] 


5 Fractional Cone Splines and Hex Splines 


We can generalize the cone splines and hex splines from Sections [3] and [4] by replacing the integer “weights” 
that represent how many times we repeat the knots in X with positive real-valued weights. 


5.1 Fractional cone splines 

We start with the fractional cone spline. 
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Figure 6 : The hex splines S( 2 ,i,i)( x | X ) and -B( 1 , 2 .i)( x | X ). The knots are x 1 = , x 2 = 

5 and — x 3 = — x 1 — x 2 = (—1, 0) T . 


Definition 5.1. Let X = {x 1 ,..., x ra ) c R s \ {0}, n ^ s, with span(X) = R s and assume a = (ai,..., a n ) e 
R" := {xeR"|a;fc > 0}. We define the fractional cone spline in terms of its Fourier transform as the 
tempered distribution given by 

__ n / | \ a k 

C a {u | X) = FT - -^ + a k n6 «w, x fe » (28) 

where 

_ f 0, a k ^ Z + 

afc “ 1 1, e Z + ' 


Our first proposal gives a closed form for C a (x | X) when n = s. 

Proposition 5.1. Let X = {x 1 ,..., x s } be a linearly independent set of vectors in R s with M = [x 1 • • • x s ], 
a = (ai,..., a s ) e R+, and {x 1 ,...,x s ) the dual basis to X. Then 


Ca( x | X) 


1 

|det(M)| 


n 


k =1 


TN 


(29) 


Proof. Consider the canonical basis E = {e 1 ,..., e s } for R s . Using (28), we see that 

» / 1 \ a k 

C a (u> | E) = FT ( --h a k 7rS (w fe ) ) 

fc=i w > 


Thus 


C'c(xlE) = 


nJL 


-b a k Tr6(u]) ) e zulXk dw. 

iu> 


When a k = 1, a k e Z and the corresponding factor is simply an a/j-fold convolution ( x k )° + * • • • (a ;*,) 3 = 
(^fc)^ -1 - When a k = 0, a k $ Z and the factor is the inverse Fourier transform of ^ ib ^ ak which is again 


13 








(30) 


P^~~y (a;fe)“ fc . Thus we see that 


C,(x E) - 1 


(Xk) 


ak -1 


(31) 


L\ r K) ' 

Note that _ _ 

CJu | X) = C a (M t u | E) . 

We use the fact that the inverse Fourier transform for f(Au:), A nonsingular, is f(Ax) in conjunction 

with (30) and (31) to infer that 

1 s / ~ k\otk~ 1 

C Q (x|X) = |det(M- 1 )|C Q (M- 1 x|E) = ra ^n^ ’ + 


|det(M)|l = l r(a fc ) 


□ 


Figure[7]illustrates the results of Proposition 5.1 ForC a (x|X) given by (291, note that supp(C a ) = cone(X) 



Figure 7: The splines C a (x|X) for X = {(1,1) T , (1, —1) T } and, left to right, a = (|, |), a = (|,2), and 

“ = (I’ 1 )- 

and it is straightforward to show that these C a (x | X) obey a recurrence formula similar to ([7]). 

Proposition 5.2. Let X = {x 1 ,...^ 5 } be a linearly independent set of vectors in and e k the kth 
canonical vector in R s . If a — e k e Ki and |a| = a\ + ■ ■ ■ + a s A s, then 

C a (x | X) = ^ <x, x k )C a _ e k (x | X). 

| Q | - a ^ 


Proof. For k = 1,..., s, 


<x, x k )C a _ e k (x I X) = <x, x k ) 


n 


c,xl> 


OLj—1 


,=i r (“j) 


/ 


<x,x fc >y ~ 2 

r (o fe -1) 


n 


<x,x J >' 


Ot-j - 1 ' 


vM r ^) J 

= (at — 1) C a (x I X). 


(ak - 1) 
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Summing both sides for k = 1,..., s, gives the desired result. 


□ 


We can generalize Proposition 3.2 for fractional cone splines on 3-directional meshes. 

Proposition 5.3. Let X be a 3-directional mesh in M 2 with x 3 = x 1 + x 2 = (1, 0) T and a e R+. Then 


C a (x | X) = ^ < 


r(ai)r(a 2 +Q 3 ) V x\ 


OL\ — 1 


1 \ Q!2 + Q'3 —1 


xi - ZjX 2 

2 / + 


_ _ X 2 

r(o:2)r(Q:i +a 3 ) ^ x% 

00 


a 2 — l 


r(ai)r(a 2 ) l S fc!r(a 3 -fc)(ai + a 2 +fc-l) I 

\k =0 


[Xi - 

(-l) fc 


) 

ir 12 ) 


2 F\ ^ 1 —ai,a2;Q2+a3;-- 

oi+a 3 -1 


2^1 l-a 2 ,ai;ai+a 3 ; 


00 


ai+ct2 + a 3 —2 


Proof. This proof follows in much the same was as that of Proposition 3.2 We can write 


C a (w | X) = C (a „ 2) (w| {x\x 2 }) • ( — + tt<j 3 6 (wi) ) , 

' VjJ\ ' 


where 0-3 = 1 if 0:3 e Z, 03 = 0 if 03 f Z, so that 

, , (®i)£ 3_1 

C Q (x | X) = <7 (ai , a2) (« | X) * ^ + 8 (x 2 ) 

-1 


x 2 > 0 
x 2 < 0 
x 2 = 0 


= f|det(M)| ]^[r(a fc ) N j f <(u, z 2 ), x 1 )? 1 0(u, z 2 ), x 2 >“ 2 1 (xi - u )“ 3 1 du. (32) 

V fe=1 / 

The dual basis vectors are given in © so that ( [32] ) becomes 

C a (x|X) = ^ ]^[ r (a fc )^ J ^u- x 2 ) Oi - u )“ 3_1 du. (33) 


If x 2 > 0, we have ^4-0 2 < ^f x 2 < x± and (|33|) becomes 


C„(x|X) = ^2 r (a fe )^ J ( U ~ x ^) ( xi-u ) a3 1 


du 


and we can use (141 to complete the proof. The case where x 2 < 0 is similar. For x 2 = 0, we use (331 to 
write 


C a ((*!, 0) | x) = U 3 n r K)j J o 1 w Q1+ “ 2 - 2 (an - u )“ 3-1 


du. 


The generalized binomial theorem can be used to complete the proof. 

Remark 5.1. If each of ai,a 2 ,oi 3 are integers, Proposition\5.3\ reduces to Proposition\3.3\ 


□ 


Figure [ 8 ] illustrates the results of Proposition |5.3| 
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Figure 8: The splines C„(x|X) for X = j (|, l) T , (— — l) T , (1,0) T | and, left to right, a = (§,§,§), 

“ = (hot’ 1 ). and « = (f>f> 2 )- 


5.2 Fractional hex splines 

In this subsection, we define a bivariate fractional hex spline implicitly in terms of its Fourier transform. We 


then use this definition to derive an explicit formula for the spline. Motivated by (21) and Definition 5.1 
we make the following definition: 

Definition 5.2. Let X = •[x 1 ,x 2 ,x 3 } be a 3-directional mesh in R 2 with x x ,x 2 linearly independent and 
x 3 = x 1 + x 2 = (1, 0) T and assume a e R+. We define the fractional hex spline of order a in terms of its 
Fourier transform (as the tempered distribution): 


B a (u | X") = (l - (l - e - l <“- x2 >)“ 2 (V<“' x3 > - l)“ 3 C a {u | X) 

where X~ = {x^x 2 ,— x 3 }. 


Taking inverse Fourier transforms of (34) gives 


B a (oJ | X') = /(x) * Cq,(x | X) 


where 


/(w) = (l - (l - e- i<u, ’ x2 >y 2 (e i<a, ’ x3 > - l)“ 3 . 

Letting z k = e k = 1,2, we use the generalized binomial theorem to write 


k i^O 


/m = E E E (I; 


= E E E 


ai \ / ai 

k\ + ty \^2 T 


{ _ 1)kl+k+ e z - klz -k 2 ( 2 iZ 2 )o 5 


= E E 

keZ 2 \ £^max(fei ,k 2 ,0) 

= E c fcllfe2 e- l <“’ M ( k -“ 3 » 

keZ 2 


£*l \ / £*2 \ [ a 3 

t-kl) V-k'2 


(_ 1 )fc 1 + fc+t \ z ~ k i z -k 2 ( ZlZ2 )°* 


(34) 


(35) 


(36) 
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where k = (fci,fc 2 ) T j ca 3 = {a 3 ,a 3 ) T , M = [x 1 x 2 ] and 


c ki,k 2 X 

^max(fei ,&2,0) 


Ot 1 


— k\) — h 


a 2 


«3 


(- 1 ) 


k!+k 2 +i 


(37) 


We have used (271 to simplify the last identity (36) above. 


We take inverse Fourier transforms of 


to obtain 


/( x ) = X c ki,kJ (x - .1/ (k a 3 )). 

keZ 2 


The relation (35) and the above argument gives rise to the following time-domain representation of B a (x X ). 


Proposition 5.4. For the fractional hex spline in Definition (5.2), we have the following explicit represen¬ 
tation: 

B a (x | X") = ^ c klM C a (x - M (k - a 3 ) | X) (38) 

keZ 2 


where the coefficients Ck lt k 2 are given by (37) and C a (x | X) is given by the relation in Proposition 
Proposition 5A is used to evaluate the hex spline plotted in Figure [9] 


5.3 



Remark 5.2. Notice that (34) can also be written in the form 


B a (uo\X-) = V£C a (u;|X), 
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where is now the Fourier transform of the fractional forward difference operator of order a := (ai, < 22 , 0 : 3 ) £ 
R+ defined by 

Vx :=nV“J, 

i= 1 

where, for any function g : R 3 —> R, 

(with a slight modification for i = 3). 


6 Hex Spline Properties and Splines of Complex Order 


We conclude the paper by stating several properties obeyed by hex splines and developing definitions of cone 
and hex splines of complex order. Unlike polynomial box splines, fractional hex splines are not compactly 
supported for a. f Z 3 . Thus we give decay rates for B a (x | X - ). We also show that hex splines are refinable 
and give an explicit formula for the refinement mask. Finally we show that a hex spline and its translates 
along the integer lattice Z 2 form a Riesz sequence and thus a R.iesz basis for an appropriate shift-invariant 
subspace of L 2 (R 2 ). 


6.1 Decay and Refinement 


We first investigate the decay rate of the fractional hex spline B a (» j X ). To this end, recall that a function 
g : R 2 —» R belongs to the Sobolev space VF r,p (R 2 ), r > 0, 1 < p < 00 , iff 


(1 


1 2 ) r/2 ff£ 


(39) 


It follows from Definition 
only on X. Hence, B a (» 


5.2 that | B a (u> 
X-) £ L 2 




K 


, where K denotes a positive constant that depends 


if |a| > 1, and B a (» I X JeL 1 


n L 2 (R 2 ) if | a | > 2. Equation 


(39) together with the above estimate on \B a (» I X )| imply that 


B a (» | X") e IW’ P (R 2 ), for r < |a| 


2 

p' 


At this point, we recall the Sobolev Embedding Theorem in R 2 [TJ Theorem 4.12]: If j e No, 1 =% m £ N, and 
mp > 2 , then 

MAJ+m^R 2 ) ^ cg(R 2 ), 

where C^R 2 ) denoted the space of all bounded continuous functions with derivatives up to order j. 

Choose m = 2 and p = 2. Then B a (» | X - ) £ Cj (R 2 ) for |a| > j + 3. In other words, for every multi-index 
p, with |/u| ^ j, | X - ) is bounded and continuous. But then [ X - ) is uniformly continuous 

and vanishes at infinity. Hence, 


B a (x | X~) e 0(x -M ), for |a| > j + 3. 
The fractional hex splines B a (» | X~) satisfy two-scale relations of the form 


B a (x | X - ) = J] h(k)B a (2x — M(k — a 3 ) | X - ), 

keZ 2 


(40) 
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with M = [x 1 x 2 ] and 0:3 = (0:3, a^) T , which are valid for |a| > 1 and for almost all x e R 2 . 

To show the validity of the above refinement equation, recall that - if such an equation holds - the coefficients 
/i(k) are the Fourier coefficients of the frequency response of the refinement filter H a . In other words, 


H a (u) = 


2B a (2uj\X~) 

B a (u>\X~) 

1 




2 u v *< 2w - xj > 


1 


,—2 


2 h 


^n(* 


3 =1 




1 — e -i<w,xi> 


(with modification for j = 3) 


1 + e 


K u 


3 \\ “3 


, (almost everywhere) 


with x 3 = x 1 +x 2 = (1, 0) T . Expanding the threefold product in the last equation above using the generalized 
binomial theorem and proceeding as in the derivation of (36), we obtain 


H = bR^T Z a ( k )' 

keZ 2 


k—a 3 )) 


(41) 


where 


«( k ) := Z 

f^max(fei ,&2 ,0) 


Oil 


l — kij \£ — k 


a. 2 


Equation (41) now yields the desired two-scale relation (40). 


6.2 Riesz basis property 


We next consider bivariate hex splines and their translates along the lattice Z 2 and show that this set 
of functions forms a Riesz sequence and thus a Riesz basis for the appropriate shift-invariant subspace of 
L 2 (R 2 ). To this end, we introduce the following fractional hex spline space. 

For the knot set X - = {x^yx 2 , —x 3 } c R 2 , where x 3 = x 1 + x 2 , and a multi-index a e R+, let 


5“ :=5“(B»(.|X-)) := /ei 2 (R 2 ) 


3a 6 f 2 (Z 2 ), /(i) = ^ a(k)5 a (x - k | X") 1. 

keZ 2 J 


(42) 


Note that for a > 1, the space S a is a principal shift-invariant subspace of L 2 (R 2 ). Our goal is to derive 
conditions on the knots such that the family {B a (» — k | X~) | k e Z 2 } forms a Riesz basis for S a . 

We proceed as follows. For B a (» | X - ), we consider the (27r, 27r)-periodic function [B a ](» | X - ) defined by 

[B a ](«|X-):= 2l^(^ + 2vrk|X-)| 2 . (43) 

keZ 2 


We employ jTf5] Proposition 5.7 (i)] to show that for |a| > 1 the family {B a (» — k | X ) | k e Z 2 } is a Riesz 
sequence in L 2 (R 2 ), i.e., that there exist two positive constants c and C such that 


0<c«S [R q ](w|X“) ^C< 00 

for almost all ueR 2 . 

For this purpose, we recall the definition of B a (uj | X~), namely, 

1 — e _l < w ’ x > \ 


( 44 ) 


B a (u X") = 


i<w, xl ) 


1 _ e -i<w,x“> 

*(o;,x 2 ) 


; i<oj,x 3 > _ y 

i(w,x 3 ) 


a G 
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Therefore, 


|5 c *(u; + 27rk|X 


-M2 


3 

n 

3 = 1 


sine 


(uj + 27rk, x- 7 )* 


(45) 


Notice that the function [S Q ,](»|X ) is symmetric about uj = so that, together with its (27t, 27t)~ 

periodicity, is suffices to establish (44) for almost all us e [0, 7 r] 2 . 


To find a positive lower bound, we consider the summand k = (0,0) in (431. As the sine function is symmetric 
about (uj, X' 5 ) = 0 it suffices to consider only positive arguments. Now, 

x J ) ^(Ixjl + lx^l) 

2 2 

If the direction x J is chosen in such a way that there exists a positive constant so that 

\x{\ + |4| 4 < 2, 

then, as the sine is decreasing over 0 < <M,xA < 7 r, 


(46) 


sine 


<M,x-?> 




—, Vw e [0 ,tt] 2 . 


Hence, we may choose c := "T[ ( —— 

• = i 

fulfilled and one may choose $3 = 1 . 


> 0. Note that since x 3 = (1,0) T , condition (46) is automatically 


To obtain an upper bound, we proceed as follows. Suppose that in addition to (46), the directions x J = 
( 4 , 4 ) T satisfy 


t , := -j $ ' 


J = 1,2, 


(47) 


where we assumed without loss of generality that 4^0- We remark on the case j = 3 below. 


Note that since uj is bounded and x 7 satisfies conditions ( |46| and ( |47| ), there exists a k 0 := (Aj 0 , 1 , 4 , 2 ) e Z 2 
so that | (uj + 27rk, x - 7 ) | = | (uj, x 7 ) + 2n(\i, x 7 )>| > 0, for all k = (fci, fo) with |fej| ^ |feo,i|» *=1,2, and all 
j = 1, 2. Let us denote the collection of all such k by K. For a k e K, the function uj >-» | (uj + 27rk, x J A | 
attains then its unique minimum value, which is strictly positive, at one of the endpoints of the square 
[0, 7 r] 2 . Denote this minimum value by 7 r(m, x- 7 )). The collection of all such m = m(k) constitutes a subset 
M of Z 2 . 


Thus 


[B a ]( W |x-)= ^ n 

keZ 2 j =1 


sine 


(uj + 27rk, 


2a-i 


= z n 

keZ 2 YJO = l 


sine 


(uj + 27rk, x^y 


2a 7 


< 


card(Z 2 \AT) + ns 

j= 1k eK 


+ zn 

keK j=l 

/ \ 2ck 

(uj + 27rk. x j y 

sine --- 

2 


[uj 


sine 


+ 27rk, x^y 


2a j 
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Now, 


(u> + 27rk, x J \ 

< 

2 



and, therefore, for j = 1 , 2 , 

z 


k eK 


(u + 27rk, 
sine —---— 


2 a? 


< 


?r| <m, x-?> | ’ 
2 

7r(2|fci| - 1) 

z 


3 = 1 , 2 ; 

> 3 = 3 > 


Vue [0, 7 r] 2 , Vke AT, 


Rewriting the last sum, produces 

1 


m eM 


z 


7r|<m,xV)| 




2 


m#(0,0) 


| (m, x-V) | 2a i 


, Jo,o) 1 |2 “' ( miin ^( 0 , 0 ) l x il 2 “ J K + r V m 2| 2Qj 


— C(Qj\oij), otj > l. 


Here, ('(Qjjoy) denotes the Epstein zeta function for the positive definite quadratic form Qj(mi,m 2 ) := 
(s{) 2 |?ni + The Epstein zeta function converges for ay > 1 pOj . 


In the case j = 3, we obtain instead 


z 

k eK 


vW 


sine - 


+ 27rk, x 3 ) 


2Q!3 


«.'2 


20-3 


2 




/ciEZ\{0} 


(2|^|-1) 2 “3 W 


z 


= C(2a 3 ), 

where £( 20 : 3 ) is the Riemann zeta function, which converges for 2a 3 > 1. 


Thus, we may select for the positive constant C as an upper bound in (44) the value 

/ 2 \ 

C := card(Z 2 \i0 + ( - J C(Qi! «i)C(Q 2; a 2 )C( 2 « 3 ) < °°i 

with ctj > 1, j = 1, 2, and a 3 > ^. 

We summarize the above findings in the next theorem. 

Theorem 6.1. Let X~ = {x^x 2 ,— x 3 } <z R 2 , where x 3 = x 1 + x 2 , he a knot set in R 2 . Suppose that the 
knots x 1 and x 2 satisfy the following two conditions: 

(i) \x{\ + \x J 2 \ <2, j = 1,2. 

xj 

(ii) -4 f Q, j = 1,2. (Assuming without loss of generality that x\ A 0.) 
x\ 

Further assume that |a| > 1. Then the family {B a {» — k [ X - ) | k e Z 2 } of bivariate hex splines is a Riesz 
basis for the principal shift-invariant subspace S a cz L 2 (R 2 ). 


It was our hope to show that in connection with the refinement equation (401 and under the hypotheses of 
Theorem 6.1 that the spaces Vf* defined by 

-L 2 (R 2 ) 


V? :=span{R a (2^.-M(k-o: 3 )|X-)|keZ 2 } ‘ ', £ e Z, 

form a dyadic multiresolution analysis of L 2 (R 2 ). The properties (a) |^J Vf* is dense in L 2 (R 2 ) and (b) 

feZ 

P'l Vf* = 0 follow immediately from [T4l Theorems 2.3.2 and 2.3.4]. Arguments similar to those employed 

.fez 

above lead to a lower Riesz bound c > 0, but it is not possible to attain a finite upper Riesz bound C. 
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6.3 Extension to Complex Orders 

Definitions |5.1| and |5.2| can be extended to complex orders z. As the previous section shows, fractional cone 
and hex splines provide continuous (with respect to order of smoothness) families of functions. A complex 
order generates a family of complex-valued functions which also contain phase information as will be seen 
below. 

We restate Definitions |5.1| and |5.2| for complex orders. 

Definition 6.1. Let X = {x 1 ,..., x™} c {0}, n 13= s, with span(X) = R s and assume z := (zi,. .., z n ) G 
C+ := {z G C" | Rezfc > 0}. We define the cone spline of complex order z in terms of its Fourier transform 
as the tempered distribution given by 


C z (u> | X) = Y\ 


= i \K“,x k ) 


1 \ z k 

+ a k nS «u;,x fc » 


(48) 


where 


® k — 


0, Zk f Z + 
1, Zk e Z + 


Similarly, we define a hex spline of complex order. 

Definition 6.2. Let X = {x^x^x 3 } be a 3-directional mesh in M 2 with x 3 ,x 2 linearly independent and 
x 3 = x 1 + x 2 = (1,0) T and assume z G C(j_. We define the hex spline of complex order z in terms of its 
Fourier transform as follows: 


B, 


; (fa> | X") := (l - e-^^y 1 _ e -i<^,x 2 >y 2 | X ) 


' ■/ 1 \ \ Z 1 / ■/ 2 \ \ z 2 / •/ 3 \ \ z 3 

1 _ e -*< w > x > \ / 1 — e _2 < w,x > \ / > — 1 \ 


ifco^ 1 ) 


i(u>,x 2 ) 


i(co, x 3 ) 


(49) 


where X = {x 1 ,x 2 , —x 3 }. 


Recall the fact that for complex numbers z and w , ur := e zlogw = e z ( lo sH+ Iar g w) j g un iq ue ]y defined for 
argw G [—7r, 7r). Hence, we consider only the principal branch of the complex exponential function. As usual, 
we set CP := 0 and w° = 1. 


Define 


and 


Ok := 0Mx fc ) := — 


1 


i(a>, x fc > 


+ a k nS «u;,x fc » 


flk := I x k ) := 


1 _ g-i(u>,x k ) 


(obvious modification for k = 3). 


i(u,x k ) 

Note that is well-defined as it never crosses the negative real axis. For Im ilk = — ^ ' xk y — = 0 iff 
(u),x k y G 27tZ and Ref^ = ^ equals 1, for (uj,x k y = 0, and 0 for ^u>,x fc )> G 27 tZ\{0}. 


With these observations, Equations (481 and (49) can be rewritten in the following form. 


C z (w | X) = C Rez (u | X) exp -^Im 2 fc - arg0 fc exp i ' l°g |Ofc| 


(50) 


22 












and 


B z (uj | X ) = -Bftez(w | X ) exp ( - Im 


Zfc • arg 


Ofcj exp U 'Yjlm.Zk ■ log|fifeM , 


(51) 


where we set Re z : = (Re zi ,..., Re z n ). 

Equations (501 and show that both the cone spline and hex spline consists of a cone spline, respectively, 
hex spline of fractional order Re z, multiplied by a modulation factor exp —^ Im Zk • arg 0*, ) and a phase 


factor exp Im Zk • log |fife| ]. A complex cone spline is plotted in Figure 


10 



Figure 10: The complex cone spline C z (u> | X) with knot set X 
(2.25 + 2.75 i, 1.5 + 3.5i). Plotted left to right are ReC z (u; | X), ImC z (u | X), and | C z (uj | X) | . 

The existence of these two factors may allow the extraction of additional information from sampled data 
and the manipulation of images. In fact, the spectrum of a complex hex spline consists of the spectrum of a 
fractional hex spline combined with a modulating and a damping factor. The presence of an imaginary part 
causes the frequency components along the negative and positive direction of x fe to be enhanced with different 
signs. This has the effect of shifting the frequency spectrum towards the negative or positive frequency side 
along x fc , depending on the sign of the imaginary part. 

Moreover, in contrast to, for instance, complex wavelet bases, the phase information exp j Imz*, • log |fife| 

V k 

is already built in, and an adjustable smoothness parameter, namely Rez, provides a continuous family of 
functions. The importance of complex-valued transforms in image analysis is discussed in [5]. 

The potential applicability of hex splines to image analysis based on the above observations will be investi¬ 
gated in a forthcoming paper and published elsewhere. 

The results in the previous subsections regarding the time domain representation of cone and hex splines, 
the decay rate and refinement property as well as the Riesz basis property easily extend to complex orders; 
a needs to be replaced by Re z and a*, by Re Zk ■ 


= {(M) 7 \( 1 >- 1 ) T } and z = 
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